シングルセルのトランスクリプトームシークエンスは、健康な状態や疾患のある状態におけるトランスクリプ トームの不均質性プロファイリングに利用できる強力なツールです。これまではショートリードシークエンス 技術が用いられていましたが、同じトランスクリプト内の距離の離れたスプライシング現象を結びつける上で 原理的に限界があるため、トランスクリプトのレベルではほとんど情報を得ることができませんでした。 Oxford Nanopore Technologies 社の技術をはじめ、ロングリードに対応したシークエンスの登場により、 シングルセルシークエンスは新時代に突入しており、トランスクリプトをアイソフォームレベルで調べることが できるようになりました。
ロングリードを用いれば、選択的スプライシングと遺伝性突然変異 のペアを簡便に検討できます。このことは、New York Genome Centre を拠点とする Dan Landau のチームが明らかにしています。 同チームは、スプライシングファクター遺伝子SF3B1 の体細胞変異 (骨髄性白血病で特に顕著に認められるドライバーの一部)とそれが トランスクリプトームに及ぼす影響との関係を調べています。スプラ イシング異常の結果、造血細胞分化障害では SF3B1 変異が発現 しますが、その正確な機序は解明されていませんでした。骨髄細胞 系列でスプライシング異常を引き起こす変異を調べることが困難で あるのは、分化の際に幅広い種類の前駆細胞が存在するだけでなく、 野生型細胞とSF3B1 変異細胞には識別できる細胞表面マーカーが 存在しないという理由があります。また、3' 末端側または 5’末端側 にバイアスが生じるショートリードシークエンスによるシングルセル 解析には、アイソフォームレベルでの変化を評価できるだけの精度 がありません
図 1 GoT-Splice ワークフローを示した図(引用元:Gaiti et al. 1 )。ナノポア 技術を用いた完全長cDNAシークエンスをGenotyping of Transcriptomes (GoT)と組み合わせることにより、シングルセルレベルで体細胞変異と 選択的スプライシングを同時にプロファイリングできます。
そこで同チームは、遺伝子発現、細胞表面タンパク質マーカー、 体細胞変異ジェノタイピングを同時にプロファイリングするための シングルセルワークフローである Genotyping of TranscriptomesSplice(GoT-Splice、図 1)を開発しました。さらに、PromethIONTM でロングリードを作成することにより、同じシングルセル内のトランス クリプトアイソフォームを調べることに成功しました1 。重要な点として、 完全長ナノポアシークエンスを用いた場合には、検出できた細胞 当たりのスプライス部位数が 4 倍になり、「ショートリードシークエンス を用いた場合のカバレッジは 3’末端側にバイアスが生じたのに対して、 トランスクリプト全体を通じてカバレッジの均一性が高く」なりました (図 2)
図 2 ERGIC3 遺伝子のシークエンスカバレッジプロット(引用元:Gaiti et al. 1 )。 ナノポアのロングリードを用いた GoT-Spliceではショートリードシーク エンス(オレンジ)と比較して、検出できる細胞当たりのスプライス部位 数が大幅に増加し(右上)(青)、スプライス部位全体を通じてシークエ ンスカバレッジの均一性が高くなります
同チームは、赤血球前駆細胞に方向づけられた細胞ではSF3B1 変異 が豊富に存在することを確認しており、SF3B1 による赤血球形成異常 の表現型に関する文献の記載と一致しています。ナノポアシーク エンスを GoT に組み込むことにより、「SF3B1 変異が細胞種に特異的 なミススプライシングを引き起こすことが示され」ました。赤血球 系列では、SF3B1 変異細胞に細胞種特異的な隠れた 3’スプライス 部位の明らかな利用が認められ、細胞周期に関連する遺伝子に 影響が及ぶことが確認されたということです 1 。
"ショートリードシークエンスは、 完全長 RNA アイソフォームとスプライ シング異常のマッピング能に限界が ある 1"
地球の反対側でも、クイーンズランド大学の Rachel Thijssen 博士が ナノポアシークエンスを用いて、血液がん、特に慢性リンパ性白血病 (CLL)と、再発時にベネトクラクスに対する抵抗性が生じる背景に ついて研究を進めています 2 。ベネトクラクスの薬理作用は、生存 促進性タンパク質 BCL2 を阻害し、最終的に白血球細胞のアポトー シスを引き起こすというものです。寛解率は高いものの、最終的には 有効性の低下が認められます。Thijssen らは、腫瘍の不均質性に ついて正確な見取り図を描くため、ベネトクラクスの単剤療法を実施 した CLL 患者から採取したサンプルをシングルセルシークエンス法に より調査しました。Landau のチームと同じく、Thijssen 博士もナノポア のロングリードを使用して変異を評価し、シングルセルのレベルで トランスクリプトの変化との関連を検討しています。その結果、アポ トーシス促進性遺伝子NOXA の非機能的トランスクリプトなど、ベネ トクラクスに対する抵抗性表現型の原因となる可能性があるものの、 これまで「正しく評価されていなかった」数々の異常なスプライシング 現象が確認されました。Thijssen の研究結果から、完全な抵抗性が 発現するのを防ぐために、ベネトクラクスは期間を限定して投与する 必要があるという考え方が裏付けられています 2。
両チームとも、体細胞変異をシングルセルのレベルでトランスクリプト の変化と結びつけることがいかに効果的であるか、完全長トランスク リプトアイソフォームを取得する上でロングリードがいかに重要である かを示しています。
"[Oxford Nanopore 社の]ロングリード シークエンスにより、同じ mRNAトランスク リプトの全体を通じてさまざまなスプライシ ング現象を定量化することも可能となり… ロングリードシークエンス独自の利点が 示された"
Single-cell transcriptome sequencing is a powerful tool for profiling transcriptomic heterogeneity in health and disease. However, in the past it has relied on short-read sequencing technologies, which provide little information at the transcript level due to the inherent limitations in linking distal splicing events within the same transcript. With the advent of long-read capable sequencing technologies, such as that provided by Oxford Nanopore Technologies, we are ushering in a new era of single-cell sequencing, where transcripts can be viewed at the isoform level.
Long nanopore sequencing reads facilitate the pairing of alternative splicing with genetic mutations. This was demonstrated by Dan Landau's team, based at the New York Genome Centre, who are exploring the relationship between somatic mutations in the splicing factor gene SF3B1 — which are some of the most prominent drivers in myeloid leukaemias — and their effects on the transcriptome. As a result of aberrant splicing, SF3B1 mutations manifest in disrupted haematopoietic differentiation, but exactly how this occurred was unknown. The difficulties in studying splice-altering mutations in the myeloid lineage are not least because of the array of progenitor cell types that are present during differentiation, but also because wild-type cells and SF3B1 mutant cells lack any distinguishing cell-surface markers. Also, single-cell analysis with 3’- or 5’-biased short-read sequencing lacks resolution in assessing changes at the isoform level.
To that end, the team developed their single-cell workflow, Genotyping of Transcriptomes-Splice (GoT-Splice), for the simultaneous profiling of gene expression, cell-surface protein markers, somatic mutation genotyping, and, due to the long nanopore sequencing reads generated on PromethION™, they could look at the transcript isoforms within the same single cell1 . Importantly, they observed a four-fold increase in the number of splice junctions per cell detected using full-length nanopore sequencing, which further ‘afforded greater coverage uniformity across the entire transcript, compared to 3’-biased coverage in short-read sequencing’.
They found that SF3B1 mutations were enriched in cells committed towards erythroid progenitors, which is consistent with the literature on SF3B1-driven dyserythropoiesis phenotype. The integration of GoT with nanopore sequencing ‘showed that SF3B1 mutations exert cell-type specific mis-splicing’. That is, in the erythroid lineage, there was obvious cell-type-specific 3’ cryptic splice site usage in SF3B1 mutant cells, affecting genes related to the cell cycle1.
On the other side of the globe, Dr Rachel Thijssen from the University of Queensland has also been using nanopore sequencing to advance her research into blood cancers, specifically chronic lymphocytic leukaemia (CLL), and the reasons underlying ventoclax drug resistance during relapse2. The pharmacologic action of ventoclax is the inhibition of the pro-survival protein BCL2, culminating in apoptosis of the leukaemic cells. Whilst this delivers high remission rates, there is an eventual loss of efficacy. Thijssen et al. investigated CLL samples treated with venetoclax monotherapy using a single-cell sequencing approach, to provide an accurate picture of tumour heterogeneity. Like Landau’s team, Dr Thijssen used long nanopore reads to assess and link mutations to changes in transcriptional read-outs at the single-cell level. They found an array of previously ‘unappreciated’ disrupted splicing events that could be attributed to the ventoclax resistance phenotype, including a non-functional transcript of the pro-apoptotic NOXA gene. Her findings support the notion that ventoclax therapy should be administered in a time-limited period to prevent the development of full-blown resistance2.
Both teams have demonstrated the power of linking somatic mutations to transcriptional changes at the single-cell level, and how long nanopore sequencing reads were pivotal in obtaining full-length transcript isoforms.
参考文献
Gaiti, F. et al. Single-cell multi-omics defines the cell-type specific impact of splicing aberrations in human hematopoietic clonal outgrowths. bioRxiv. 495292 (2022).
Thijssen, R. et al. Single-cell multiomics reveal the scale of multi-layered adaptations enabling CLL relapse during venetoclax therapy. Blood. 2022016040 (2022).